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The ekpyrotic and cyclic universe scenarios have revived the idea that the density perturbations 
apparent in today's universe could have been generated in a 'pre-singularity' epoch before the big 
bang. These scenarios provide explicit mechanisms whereby a scale invariant spectrum of adiabatic 
perturbations may be generated without the need for cosmic inflation, albeit in a phase preceding 
the hot big bang singularity. A key question they face is whether there exists a unique prescription 
for following perturbations through the bounce, an issue which is not yet definitively settled. This 
goal of this paper is more modest, namely to study a bouncing Universe model in which neither 
General Relativity nor the Weak Energy Condition is violated. We show that a perturbation which 
is pure growing mode before the bounce does not match to a pure decaying mode perturbation 
after the bounce. Analytical estimates of when the comoving curvature perturbation varies around 
the bounce are given. It is found that in general it is necessary to evaluate the evolution of the 
perturbation through the bounce in detail rather than using matching conditions. 
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, In the inflationary Universe scenario, large scale structure is generated by the stretching of subatomic scale quantum 
' fluctuations to macroscopic length scales, during an epoch of superluminal expansion supposed to have occurred prior 
to the radiation dominated era. Indeed such expansion appears to be essential to the creation of correlations on 
super-Hubble radius scales, a feat prohibited by causality in the standard hot big bang era. 

However, cosmic inflation does not resolve the problem of the initial singularity, and it seems clear that a more 
complete theory than inflation is needed to deal with it. But if time existed before the the initial singularity, it is a 
logical possibility that the large scale structure and density perturbations we see today were generated during this 
\Q , pre-singularity epoch. The ekpyrotic and cyclic scenarios [1,2] provide explicit realisations of this idea (pre-figurcd in 
the 'pre-big bang' scenario of Veneziano et al. [3]). In the ekpyrotic and cyclic scenarios, scale invariant large scale 
density perturbations are generated during a phase which is contracting (in the Einstein frame) before the bounce 
to expansion [1]. The task of matching these perturbations across the bounce to the expanding epoch is highly 
challenging, since general relativity must break down there, and the use of string theory methods will ultimately be 
■ essential. There are indications that for the particular type of singularity involved here (namely the collapse of a single 
O ,' extra dimension) the divergences are relatively weak [4] , but the matching issue remains unsettled at the present time 
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i-C ' In a contracting universe, the growing mode adiabatic density perturbation corresponds to a shift in the time to the 
, 'big crunch'. A perturbation in the space curvature, however, is a decaying mode perturbation because it becomes 
increasingly irrelevant as the scale factor shrinks. (These statements both hold in Einstein frame and comoving gauge). 
However, in an expanding universe, the situation is reversed. A time delay is now a decaying mode perturbation and 
the curvature perturbation is the growing mode. This has led several authors to suggest that in the ekpyrotic and 
cyclic scenarios, where scale invariant density perturbations are generated in the collapsing phase growing mode, that 
these perturbations would match on to pure decaying mode perturbations in the subsequent expanding phase. Indeed 
this result is obtained if one insists on matching the curvature perturbation across the bounce [5]. 

However, this behaviour would appear surprising from a more physical viewpoint. In the models being considered, 
gravity is attractive throughout. The growing mode instability is driven by gravitational attraction and it is physically 
implausible that it should precisely reverse at a bounce so as to dissipate after the bounce and render the final universe 
perfectly homogeneous. 

In this paper we address this issue by considering a model in which the issue can be definitively settled within 
conventional general relativity. It is well known that a massive scalar field in a closed Friedman-Robertson- Walker 
(FRW) Universe can lead to a bounce without a singularity, and without violating the Weak Energy Condition. This 
model has a long history starting with Schrodinger [8] , however he did not include the back-reaction of the scalar field 
on the background scale factor. It was also employed in semi-classical studies of quantum effects (see for example [9]), 
and of course is a common model in quantum cosmology and studies of the no boundary proposal (see for example 
[10], and more recently [11]). The classical dynamics for the spatially homogeneous case have also been investigated 
quite extensively, see for example [12]. 

A preliminary study of cosmological perturbations in such a Universe (and in other bouncing models) was performed 
by Hwang and Noh [13], with inconclusive results since the usual large scale approximations break down around the 
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bounce. We develop new approximations that do not break down around the bounce, and using numerical simulations, 
accurately determine the propagation of linear perturbations through the bounce. We show that growing mode 
perturbations developed in the collapsing phase do not match to pure decaying mode perturbations in the expanding 
phase, thus there is no contradiction with the physical argument given above. Of course, this does not at all prove 
that the same is true in the ekpyrotic/cyclic models: as the original papers made clear, an unambiguous matching 
condition is still required. But the present work does show that propagation of a growing mode perturbation across 
a bounce is possible at least in this toy model. 

We would also like to briefly mention recent work on the ekpyrotic/cyclic scenarios. Durrer and Vernizzi [14] showed 
how matching across the bounce depends on what surface the matching is done on and the intrinsic stress energy 
of that surface. They showed that only for the special case of matching on constant energy density surfaces with 
no surface tension that the growing mode in the collapsing phase matches completely onto the decaying mode in 
the expanding phase. Generally, some of the growing mode in the collapsing phase will go to the growing mode in 
the expanding phase. We should also mention a recent study of the possibility of isocurvature perturbations in the 
Ekpyrotic model [15]. Peter and Pinto-Neto [16] studied the matching problem in a hydro-dynamical fluid, and in a 
second publication employed a scalar field with negative kinetic energy in order to obtain a regular bounce [17]. The 
problem of matching quantum fields across a bouncing universe of the type relevant to the ekpyrotic/cyclic scenarios 
has been considered by Tolley and Turok [18], and a number of studies have been made of analogous bouncing models 
within string theory [19]. 



II. MODEL 



The equation for the scale factor in a FRW Universe is given by 

1 

6M P 



a = 2 (p + 3p)a (1) 



where a is the scale factor, dot denotes a derivative with respect to time, p is the density, p is the pressure and 
Mpi = 1/V87TG is the reduced Planck mass. As the scale factor is always positive and d is negative for a collapsing 
Universe, it follows from Eq. (1) that the Strong Energy Condition (p + 3p ^ 0) must be violated for a > 0. Which 
is what is needed for a bounce at a finite scale factor to be possible. This is the same as the condition for inflation. 
It is satisfied by a scalar field that is potential dominated as 

P=\i 2 + V (2) 

and 

V=\^ 2 -V (3) 
where </> is the value of the scalar field and V is its potential. For simplicity we shall only consider a quadratic potential 



V = l -m 2 4> 2 (4) 



where m is the constant mass of the scalar field. 
Using the Friedman equation 



H 2 = -^p-^ (5) 
3M P1 2 a 2 y ' 

it can be seen that the background curvature (K) must be positive for the Universe to bounce. This is because at the 
bounce H = and p > 0. As p + p > 0, it follows that the Weak Energy Condition is never violated by this model. 

Although the above arguments show that it is possible for a closed Universe with a scalar field to bounce, it does 
depend on the initial conditions. The scalar field evolution is given by the Klein-Gordon equation 

4> + 3H<fi+V <p =0 (6) 

where = dV/d(f>. In a expanding Universe H > and so the 3H(j) acts as a friction term. While in a collapsing 
Universe H < and so the 3H<f) acts as an anti-friction term. 
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So, it is possible that the the system will always be kinetic energy dominated and go to <fi = ±00 and a = instead 
of bouncing. It is only for a small subset of initial conditions that such a Universe would bounce [12]. The bounce 
will be symmetric in the background quantities if <p — when a = 0. We use this case to study the evolution of 
perturbations through the bounce which is taken to occur at t — 0. Figure 1 summarises the sequence of events. 




FIG. 1. A sketch of the background dynamics of the model investigated in this paper. The scalar field is represented by the 
ball and ln(a) is proportional to the size of the ball. The scale factor changed by approximately 15 e-folds from its minimum 
to maximum size in the numerical integration. First the scalar field rolls up the potential as the Universe collapses. Then 
when the Universe reaches a non-zero scale factor ao, it starts to expand again and roll down the potential. The background 
dynamics are time-symmetric around the bounce. 

Around the bounce we can make the approximation <fi w and (f> = <fro where subscript is taken to the value of 
the quantity at the bounce (t = 0). Using Eqs. (1), (2) and (3) this gives 

a w ao cosh(i/ao). (7) 

This is the same as the equation for a closed de Sitter space. If we choose K = 1 then from Eqs. (4) and (5) 
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The Hubble parameter around the bounce is then given by 

H = — tanh(t/a ). 
a 



(8) 



(9) 



So the time scale around the bounce is a - In order to approximate the evolution of the scalar field around the bounce 
we rewrite Eq. (6) as 



(<M 3 



Integrating both sides gives 



-V> 3 . 



U^a 3 dt 



where (<^>)o = was used. Using w (U^)o and Eqs. (4) and (7) we get 

- In(cosh(t/a0)) - 1 ) + <f> - 



1 2 2 , ( 1 



cosh (t/a ) 



(10) 



(11) 



(12) 
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Away from the bounce (t > do), this equation can be approximated by 



1 + ln(2) -■!-!■ + 1 



(13) 



This is the same equation that is obtained by using the slow roll approximations to Eqs. (5) and (6) (see for example 
[20]) 



zh<i> + u « o 



and 



ff 2 



3M P1 2 



rV. 



Which will be a good approximation when the slow roll parameters [20] 



Mp, 2 /W 2 
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and 



= M PI 2 - 



d 2 U 



(14) 



(15) 



(16) 



(17) 



are small and there has been enough expansion to make the K/a 2 in Eq. (5) negligible. For the quadratic potential 
given in Eq. (4) both conditions become 



<t> 2 > 2M Pl 2 . 

For \t\ > oq, a good approximation for H is obtained by substituting Eq. (13) into Eq. (15) 



H(\t\>a ) 



sign(t) 



a 



M, 



pi 



-1- — +ln(2) + 1 



(18) 



(19) 



This can then be integrated to obtain an expression for a(t) and the integration constant can be found by matching 
to the bounce approximation Eq. (7) 



. ^ \ «o f\t\ /, lA , n \t\\ M P1 2 \t\ 
2 \ a o V a J (pQ a 

If we further assume |</>o|/Mpi » |t|/a then Eqs. (13), (19) and (20) can be simplified to 



^(1*1 > «o) ~ <^o, H(\t\>a ) 



sign(i) dH 
~~dt 



M P / 1 



a 



l*l>oo 



2 2 , o(|t| > oq) « — exp 
& a l \do 



(20) 



(21) 



In which case ao is also the time scale away from the bounce. 



III. PERTURBATION EQUATIONS 

For studying the evolution of the perturbations we use the Bardeen formulism [21,22]. For purely scalar perturba- 
tions in the zero-shear (conformal Newtonian) gauge, the perturbed metric is given by 



ds 2 = -(1 + 2$)dt 2 + a 2 (t)(l - 2^)-y ij dx i dx j 
where x % is the comoving spatial coordinate and for a closed Universe 



Jij = S i3 ( 1 + -Kx k x k 



(22) 



(23) 
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where Latin indices range from 1 to 3, $ and "J - are the metric perturbations and Sij is the Kronecker delta. For a 
scalar field there is no anisotropic stress and so there is the further simplification "J" = $ [22]. 
The equation for the evolution of $ can be derived from the Einstein equations [22,23] and is 

$ + - $ + ^(-V 2 - 4K")$ + 2 - Hjj $ = (24) 

The usual harmonic decomposition is performed on the perturbations so that in a closed Universe the eigenvalues of 
— V 2 arc (fc 2 = [n 2 — 1]K) [24] where n > 3 is an integer. The n = 1 mode is homogenous and the n = 2 mode is a 
pure gauge mode [21]. We will be looking at the evolution of an individual mode with time. So the dependence of <f> 
will be in terms of n and t. 

An approximation for $ around the bounce (|t| < do) can be found by substituting Eqs. (12) and (7) and K = 1 
into Eq. (24) and working with the harmonic components of <& 

d+ Aanhft/ ao ) _ 2 V + n 2 -5 $ 
y «o 0y a 2 , cosh (t/ao) 

+2 / J_ _ tanh 2 (£/a ) _ tanh(t/a ) «A $ = Q 
The equation can be reduced to a solvable form by a non-linear transform of the time coordinate 

"- h ( tanh G£-r)) (26) 

where z = V^I. Then Eq. (25) becomes 

, 6(-l + e*0(l + e *Q gg / n2 5 4(1 + e 2 -) 2 \ 
9l/ 2 ^ (-l + 4e 2 --e 4 -)(l + e 2 -) dv \ -1 + 4e 2v - e*» J { ' 

Which has the solution 

$(|t| £ ao) « Ci*i + C 2 *J (28) 

where Ci and C2 are integration constants, superscript * denotes the complex conjugate and 

= e (^> [_( n _ i)( n _ 2)e 6,y + 3(n + l)(n - 2)e 4,y , > 

+ 3(n + 2)(n - l)e 2y - (n + 2)(n + 1)] . [ ' 

We will use this solution to check our numerical calculations in Section IV. 

For \t\ > a the Universe is effectively flat (K = 0) and we can use the the background approximations given in 
Eqs. (13), (19) and (20). As seen from Eq. (24), we can also treat the perturbations as if we were in a flat Universe 
provided n 2 — 1 ^> 4. In which case the usual transformation can be used [22,25] 

?i = (30) 

which has the equation 

u" + (fc 2 - 6"/6)u = (31) 
where prime indicates the derivative with respect to conformal time r with adr = dt and 

e = - (32) 

a<p 

The 'large scale condition' is 

fc 2 < e"/6. (33) 
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A simple expression for the time spans that the large scale condition is satisfied can be obtained by assuming that 
\(/)o\/Mpi 3> \t\l a o an d using Eqs. (21), (32) and (33) which gives 

aO V M Pl J v ' 

As is expected the larger k the further from the bounce the large scale approximation becomes valid. Also, the large 
scale condition is always violated for \t\ > a as k 2 >= 3 2 — 1 and \<fr \ > Mp\. When the large scale condition is not 
satisfied the solutions of Eq. (30) and therefore Eq. (24) are clearly oscillatory with frequency k. 

If the large scale condition is satisfied then we can approximate the solution to Eq. (31) by a truncated series of k 2 : 

oo 

u = J2 C Mk 2 >. (35) 

1=0 

Substituting this into Eq. (31) and matching coefficients on the left and right side gives the zeroth order approximation 

u = C 3 9 + C A J 1 dr. (36) 
Using Eqs. (13), (19), (20), (30), (32) and (36) and redefining the constants C 3 and C 4 gives 

$(|t| £ ao) « C 3 ^ + d (1 - v^erfcx(VO) (37) 



where 



and 



erfcx(^) = erfc(^) (39) 

is the scaled complementary error function. If we further assume |0o|/Mpi ^ |i|/«o, then substituting Eq. (21) into 
Eq. (37) gives [26] 

*(|*|/ao > 1) « 2 fc 3 sign(t)^exp (-\t\/a ) + C 4 --^-) . (40) 

Although the form of the solution (Eq. (37)) is the same before and after the bounce, it is not necessary that the 
integration constants (C3 and C4) should be the same. The slow-roll approximation for Eq. (24) is not the same 
before and after the bounce due to the presence of the absolute value signs in Eqs. (13) and (19). So we will have 
four integration constants: 

c rc ?ow , tZ-a 

{ C docay' ^ ~ a 



and 



<\ { ( 'rr-- J-" 00 (42, 



The + and — superscript indicate after and before the bounce. The 'grow' and 'decay' subscripts indicate that the 
constant is the amplitude of the growing and decaying mode respectively. 

In principle there should be two constraint equations relating the four integration constants as Eq. (24) is second 
order and so should only have two integration constants. That is, there should be a matching rule expressing C c j ecay 



and C+ row in terms of C decay and C grow . 
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IV. NUMERICAL SOLUTIONS 



Although Eqs. (28) and (37) provide good approximations around and away from the bounce. They do not overlap 
sufficiently to provide a solution for the whole domain of interest for all initial conditions. In principle, the stage 
before the global curvature becomes important, \t\ > ao, and when the large scale condition is violated, Eq. (34), could 
also be included in the matching. But the analytical approximations are found still not to be sufficiently accurate in 
the overlap regions for general initial conditions. Therefore it is also necessary to solve Eq. (24) numerically. 

Before doing a numerical integration it is important to identify if there are any singularities in the coefficients of 
Eq. (24). A series solution for Eqs. (1) and (6) is 



and 



a = ao 



1 + 5U' + " (41 



Substituting Eqs. (43) and (44) and r = t/a into Eq. (24) gives 

<9 2 $ 4 9$ A -n 2 + 5- 

_ 



2r 2 



(2 + r 2 )r dr 



(2 



$ = 0. 



(43) 



(44) 



(45) 



As can be seen there is a singularity in the coefficient for <9$/<9t at r = 0. This means that it is not possible to 
numerically integrate Eq. (24) through the bounce. However, we can solve Eq. (45) to give an analytical solution very 
close to the bounce 



$ = {2 + t 2 )- n 
{C 5 F(-A+±(-l 



33), -A - 1(1 + V33),- 



2 ' 2 ' 1 



+ C 6 r 3 F (-N + 1(5 - V33), -A + |(5 + V33), §, -±r 2 ) } 



(46) 



where F(-) is the hyper-geometric function, A = yj (n 2 — l)/2 and C5 and Cq are integration constants. A series 
approximation can be used to evaluate Eq. (46) around the origin 



^ Ccvcn^cvcn + Codd^odd 

where C even and C dd are integration constants and 

4>cven = i(n 2 -5) (-) +1 + 0(4) 

and 



$odd = ( — J +0(5). 



(47) 



(48) 



(49) 



Eqs. (47), (48) and (49) can then be used to generate initial conditions on either side of the bounce for an adaptive step- 
length Rungc-Kutta numerical integration routine [27]. We set these initial conditions for the numerical simulation 
at t/a = ±1CP 4 . The relative accuracy for the Runge-Kutta routine was set to 1CP 14 . The numerical results were 
found to be insensitive to the precise values used. 

The parameters used for the simulations were <fio — — 50Mpi and m = 10~ 6 Mpi. As the system of equations Eqs. (1), 
(6) and (24) only depends on m through the combination mt, changing m will just change the scaling of time. As can 
be seen from Eqs. (26), (28) and (29), around the bounce the only dependence of <po is through the dependence on a , 
see Eq. (8). So around the bounce changing </> will just change the time scale of the problem. From Eqs. (21) and 
(40), it can be seen that increasing <p will effect the time scale through its effect on ao (Eq. (8)) and the normalisation 
of the C3 and C4 modes. 

The only other parameters that need to be specified are the initial values for $ and $. The numerical solution 
for these can be found by choosing C OV on and C dd in Eq. (47) so that the numerical extrapolation has the required 
values of $ and $ at the initial time. This technique is known as shooting [27]. 
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Figure 2 shows a plot of $ for the wave- number corresponding to n = 3. The initial conditions at t/ao — — 20 were 
set to be a pure growing mode $ = H/a using the shooting method. The decaying mode in the collapsing phase and 
the growing mode in the expanding phase was evaluated by matching Eq. (37) at the initial and final times. They 
are plotted for those times for which the large scale condition, Eq. (34), is satisfied. This corresponds to \t\/ao > 5. 
Due to the finite accuracy of the shooting algorithm, there is still a finite amplitude of decaying mode a t/ao = — 20. 
However, by the time the large scale condition, Eq. (34), is violated, the decaying mode is negligible relative to the 
growing mode. 

io- 9 

10 -io 

io- 11 

1Q- 12 

io- 13 
io- 14 
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io- 17 

10 -18 

10" 19 

10 -20 
10- 21 

-20 -10 10 20 

t/a 

FIG. 2. A plot of the absolute value of the zero shear curvature perturbation $ for the wave number corresponding to n = 3. 
The solid lines are for positive $ and the dotted line is for where is negative. The collapsing phase decaying mode and 
expanding phase growing mode (Eq. (37)) are plotted as dashed lines. 

As can be seen from Eqs. (47), (48) and (49), $ and its first and second time derivatives are finite at t = 0. This 
shows that there are initial conditions that lead to linear perturbations around the bounce [7,16]. 

Figures 3 and 4 show examples of numerical solutions and the bounce approximation, Eq. (28), for small and large 
wave numbers. The initial conditions for Figure 3 are the same as those for Figure 2. In Figure 4, the initial conditions 
are set so that at t/ao — ~ 20 there is a pure decaying mode, Eq. (40), $ = Mpi /<f%. It is an important confirmation 
of the numerical integrations that they match the bounce approximations so well for |£|/a ~ !• 
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FIG. 3. A plot of the absolute value of the zero shear curvature perturbation ($) and its approximation around the bounce, 
Eq. (28), for n — 3. The solid lines are for positive $ and the dotted line is for where $ is negative. The absolute value of the 
approximation is plotted as a dashed line. 
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FIG. 4. A plot of $ (solid line) and its approximation using Eq. (28) (dotted line) for n — 10. 

V. COMOVING CURVATURE 

In inflation, the use of the curvature on constant density hyper-surfaces (Cbst) [28] or the curvature on comoving 
hyper-surfaces (TZ) [29] is very useful because in an expanding flat Universe with no decaying mode or entropy 
perturbation they are constant on scales larger than the Hubble horizon (1/H). We will be examining whether this 
can be used for the present model. 

The comoving curvature perturbation for any background curvature (K) is [30] 
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ft = * + 2M m 2 H — (50) 

p + p 

where for a scalar field = Note that the formula for the comoving curvature giving in Mukhanov ct. al. [22] 
(Cmfb) only applies for the flat case (K = 0). Alternatively the curvature on constant density hyper-surfaces [28] 
could be used. However, away from the bounce, when the decaying mode is sub-dominant, they are all approximately 
equal up to a change in sign and so they will all serve the same purpose below. 
Using Eqs. (43), (44), (47), (48) and (49) in Eq. (50) gives 

TZ(0) = C cvcn (2n 2 - 7) (51) 

which shows that TZ is finite at the bounce. 

Using Eqs. (2), (3), (4), (13), (19) and (37) into Eq. (50) gives 

TZ w C 4 . (52) 

From Eq. (52) it appears that TZ = when the large scale condition, Eq. (33), holds. However, Eq. (52) was derived 
using the zeroth order truncation of fc 2 , Eqs. (35), (36) and (37). As is well known [22,25,20], TZ will not be constant 
if the series expansion, Eq. (35), is extended to order fc 2 . Alternatively, the time change of R can be found using the 
general equation [20] 

TZ = (53) 
P + P 

where <5p C omov is the pressure perturbation on comoving hyper-surfaces. For a single scalar field this becomes [31] 



■ H k _ 

(54) 



As a shrinks exponentially as the bounce is approached, Eq. (20), TZ may no longer be approximately constant. Thus, 
71 can vary on large scales as the pressure perturbation on comoving hyper-surfaces grows exponentially as the bounce 
is approached. 

The pressure perturbation can be decomposed as [20] 

P 

Sp = -Sp + Sp nad (55) 
P 

where Sp is the density perturbation and <5p na d is the gauge invariant non-adiabatic component of the pressure 
perturbation. For the quadratic potential in this paper, both the adiabatic and non-adiabatic components of the 
pressure perturbation on comoving hyper-surfaces grow exponentially as the bounce is approached. Where as for 
a massless field, such as in the stage just before the bounce in the Ekpyrotic scenario [1] there is no non-adiabatic 
component to the pressure perturbation [31] and it is only the adiabatic part of the pressure perturbation that grows 
and causes TZ to vary [7]. 

Once the large scale condition, Eq. (33), is violated then $ and hence TZ will become oscillatory. In order to quantify 
whether TZ varies before this , we look at whether there is a significant relative change in TZ in one Hubble time which 
is approximately given by the time when 



1 TZ 
HTZ 



> 1. (56) 



Assuming the C4 mode is dominant in Eq. (40) and substituting Eqs. (21), (40), (52) and (54) into Eq. (56) gives the 
condition for significant variation for TZ to be 

4fc 2 |ft|cxp(-2|t|/a ) > \TZ\ (57) 

which implies there will be significant variation in TZ due to the comoving pressure perturbation even on large scales 
for 

^ < ln(2fc). (58) 
ao 



10 



However, comparing this with the large scale condition, Eq. (34) and remembering that we have assumed \<po\ ^ Mpi, 
it is clear that the large scale condition is always violated before the condition in Eq. (58). So, we can conclude that 
if the C4 mode is dominant in Eq. (40), then the large scale condition will be violated before TZ is effected by the 
comoving pressure perturbation. 

Assuming the C3 mode is dominant in Eq. (40) and assuming |(/> |/Mpi 3> |i|/ao, substituting Eqs. (21) and (40) 
into Eq. (54) gives 



i* 



-2$*fc 2 



Afpi' 



■ exp 



«0 



a 



(59) 



where $» = $(t«) and t* is arbitrary except for satisfying the large scale condition Eq. (33). Substituting Eqs. (59) 
and (52) into Eq. (56) gives the time period for which TZ will not vary on large scales due to the comoving pressure 
perturbation 



1*1 



>- 

ao ~ 3 



In 2k 2 



Mi 



pi 



ao 



(60) 



If t* is taken to be the time when the large scale approximation breaks down, Eq. (34), then using Eqs. (40), (52) and 
(60) implies the condition for the comoving pressure perturbation to cause TZ to vary before the large scale condition, 
Eq. (34) is violated is 



]*s| > 1 
|$ 4 | 4 V Mpi 



k 



(61) 



where $3 and $4 are the C3 and C4 modes in Eq. (40). It follows that provided the $3 mode is sufficiently larger 
than the <I>4 mode, TZ can vary significantly even while the large scale condition, Eq. (34), holds. However, for large 
enough k the modes will always violate the large scale condition before the comoving pressure effects 7Z. 

Eq. (60) can be tested against the numerical results in Figures 2 and 5. Separate values of and C4 are needed 
for before and after the bounce. The values of $* can be read off for some |t*| of Figure 2 which satisfies the large 
scale condition, Eq. (34) and is where the growing mode is dominating. The values of C4 can be evaluated from the TZ 
constant parts of Figure 5. These values can be substituted into Eq. (60) to determine the region in which the entropy 
perturbation causes R to vary. These limits are drawn in with dashed vertical lines in Figure 5. As can be seen they 
accurately delimit the region in which TZ changes. From Eq. (34), it can be seen that the large scale condition holds 
for |t|/ao ^ 5 and so TZ is varying in the collapsing phase even when the large scale condition holds. 




FIG. 5. A plot of the absolute value of the comoving curvature perturbation (TZ) corresponding to the zero shear curvature 
perturbation (<£>) in Figure 2 which was for n = 3. The solid line is positive values of TZ and the dotted line is negative values. 
The vertical dashed lines enclose the time domain in which TZ is expected to vary according to Eq. (60). 
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Eq. (52) can be used to write Eq. (37) as 

S(|t| > |t t |) w C 3 ^ +7e [1 - ^Verfcx(^)] . (62) 

Where |tf| satisfies both Eqs. (34) and (60), i.e. it delimits the time on which 1Z can vary. The even and odd modes 
obtained from the numerical solution in Section IV can be decomposed into growing and decaying parts as in Eqs. (37), 
(41) and (42) 

g. _ J C- vca ^ ow H/a + fteven(-|*tl) t 1 ~ V^erfcx^)) , t < -|i t | (f ,^ 

CV ° n ~ I C+ cn , dccayJ ff/a + ft cvcn (|t t |) (1 - V^erfcx(VO) , * > |t t | ^ 

A similar equation holds for the odd mode 

* „ J Codd.grow^/ + ^odd(-|t t l) (1 - ^erfcx(^)) , t < -|t t | .... 



(65) 



Co + dd,dcca y tf/« + ^odd(|*tl) (1 " V^erfacty)) , t > |t t | 
The growing and decaying mode also have even and odd properties: 



H 

a 



H 

a 



and 

(1 — y/TTipericx(ip))\ t = (1 — y / 7n/' er f cx W'))|_ t • (66) 

The fact that the growing and decaying mode have even and odd properties does not imply they do not both contribute 
to $evcn and $odd- It does mean that it is possible to reduce the number of unknown coefficients used to describe the 
even and odd modes. From Eq. (63) 



C~ * 

^even,grow ^ 



+ n cvcn (-m\) (1 - \/W>erfcx(V0)|_ t = Cven.decay f 
-t a 



from which it follows that 



C~ * 

^even,grow ^ 



r+ — 

even, decay 



+ n 

even (|t t |) (l-^erfcx(V0)| t (67) 



(68) 



which when combined with Eq. (65) implies 

^even,grow ^even. decay' (69) 

Also from Eq. (67) it follows that 

^cven(-N) (1 - V^erfcxW>))|_ 4 = fteven(|* t l) (1 ~ V^erftX^))^ (70) 

which when combined with Eq. (66) implies 

ftcvcn(-|*tl) =^cvcn(|tf|). (71) 

Similarly, using Eqs. (64), (65) and (66) it follows that 

(j- — rj+ (jo) 

odd, grow odd, decay' v / 

and 

^odd(-|*tl) = -^odd(l*tl)- ( 73 ) 
An initially pure growing mode in the collapsing phase could be constructed as follows 

Sr ° W ^odd(-|*tl) ^even(-|i t l) 1 ' 
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provided 7£ dd(— |*tl) an d ^even( — |*fl) are non-zero. Using Eqs. (63), (64), (71) and (73), Eq. (74) becomes 



a V7£ ddH*tl) ^ ov „»(-|t t |) 



grow 



H 



Kdd(-|t t D ^ ovo „(-|t t |) 



*<-|*tl 

- 2(1 - V^erfcx^)), t > |t t | 



(75) 



Substituting Eq. (74) into Eq. (50) gives 



ft, 



■dd 



gl '° W ^odd(-|i t l) ^cvcn(-litl) 

and substituting Eq. (75) into Eq. (50) gives 



Tla 



0, t<-|*tl 
-2, * > 



(76) 



(77) 



As we are free to multiply Eq. (74) by a k dependent constant, it follows that lZ gIOW (t > \tj\) may be some other 
value besides minus two. But it is only zero if $ grow = 0, i.e. no perturbation at all. Therefore we have shown 
analytically that an ingoing pure growing mode leads to a non-zero outgoing growing mode provided both the odd 
and even modes have a C4 component. However, we can not say what the scale dependence will be without knowing 
the scale dependence <fr gr0 w Figure 6 shows a plot of the various quantities in Eq. (76), where — \t^\ is taken as the 
minimum time in the figure. 



-1 - 



-2 - 



-10 





t/ a n 



1 



FIG. 6. A plot of the comoving curvature perturbation 1Z for n = 3. The initially pure growing mode as constructed in 
Eq. (76) is the solid line. The normalised (as in Eq. (76)) even (dotted line) and odd (dashed line) are also displayed. 



This numerical integration is in agreement with Eq. (77). 



VI. CONCLUSION 



In this article we have examined the way perturbations evolve across a bounce in a scalar field dominated closed 
Universe. The background model naturally bounces without violating General Relativity or the Weak Energy Condi- 
tion. Analytical approximations were found for the background and perturbations. The perturbation approximations 
were found to work well around the bounce \t\ ^ ao and away from the bounce in the slow-roll regime \t\ £ ao- 
However the overlap between the approximations was not sufficient to use the analytical approximations alone for 
arbitrary initial conditions. But they were useful in checking and interpreting the analytical solutions. The numerical 
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solutions were complicated by the ordinary differential equation for <I> having a singular coefficient at the bounce. 
Thus a Taylor approximation had to be used around the bounce and shooting methods used to obtain the numerical 
solution for particular initial conditions for a wider domain. The Taylor expansion around the bounce could also 
be used to show that the perturbations and their first and second time derivatives are finite at the bounce and so 
non-linearities in the Einstein tensor are not inevitable, but depend on the choice of initial conditions. 

It was shown that all wave number modes violate the large scale condition for \t\ > ao and so become oscillatory 
around the bounce. We also found during the collapsing phase the growing mode could cause the comoving pressure 
to increase sufficiently to make 1Z vary significantly as the bounce is approached. The same is true for the decaying 
mode as the bounce is approached from the expanding phase side. Thus, in general we do not expect 1Z to be constant 
around the bounce. A formula predicting when 1Z starts to vary was given. We also showed numerically that, for this 
model, a pure ingoing growing mode (1Z~ = 0) did not lead to a pure outgoing decaying mode (1Z + = 0). 

This model shows that in general the dynamics of the bounce need to be taken into account when evaluating how 
perturbations change across a bounce. 
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